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ABSTRACT 


On the basis of A. T. Young’s suggestion that elemental sulfur is respon- 
sible for the yellow color of Venus, we compare calculations at 3. A ym of the 
reflectivity phase function of two sulfur containing inhomogeneous cloud 
models with that of the homogeneous model of 'Hansen and Hovenier. The inhomo- 
geneous drop model has a spherical core of solid sulfur within each sulfuric 
acid drop, while the inhomogeneous cloud model has a bimodal size distribution 
with irregular sulfur particles comprising the second size mode near 13 pm. 

The calculations were made at 3.4 pm because the ability to detect small con- 
tributions from non-absorbing species to the sunlight scattered by the Venus 
clouds is greatly enhanced (~40 times) in the infrared absorption band of 
sulfuric acid. Assuming reflectivity observations with 25% or less total 
error, comparison of the model calculations leads to a minimum detectable mass 
of sulfur equal to 7% of the mass of sulfuric acid for the inhomogeneous drop 
model. For the inhomogeneous cloud model the comparison leads to a minimum 
detectable mass of sulfur between 17% and 38% of the mass of the acid drops, 
depending upon the actual size of the large particles. That amount of sulfur 
is equivalent to detecting about one large sulfur particle per 10^ sulfuric 
acid drops. We conclude that moderately accurate 3.4 pm reflectivity observa- 
tions are capable of detecting quite small amounts of elemental sulfur at the 
top of the Venus clouds. 
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*Tel-Aviv University 
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INTRODUCTION 


It is generally well accepted that sulfuric acid is the principal constit- 
uent of the Venus clouds, and that the planet’s yellow color is caused by 
absorption of the -blue end of the reflected solar spectrum. It is also known 
that sulfuric acid is transparent in the blue and visible, and is thus incapa- 
ble of producing the yellow color of the planet. Such reasoning suggests "that 
an additional cloud constituent be considered . 

There are three ways in which an additional blue absorber could be added 
to the Venus cloud aerosol. The first has the additional constituent com- 
pletely dissolved in the liquid sulfuric acid cloud drops , and the second 
places the constituent as a solid core inside each liquid drop. In the third 
way, particles (or drops) of the blue absorber are mixed inhomogeneously among 
the acid drops of the clouds. A fourth way, which we shall not treat here, 
adds the blue absorber to the cloud as an unknown gas . 

Recently, A. T. Young (ref. 1) amplified an earlier suggestion by Hapke 
(ref. 2) that elemental sulfur is the additional absorber because it has a 
temperature sensitive blue absorption edge capable of explaining the visual 
color of Venus. If sulfur is the new constituent, it is not added to the 
clouds in the first way described above because sulfur does not dissolve in 
sulfuric acid (personal communication, 0. B. Toon, Ames Research Center, 1978). 
However the other two ways, the inhomogeneous drop and the inhomogeneous cloud j 
are distinct possibilities. Young (ref. 1) has already proposed an 
inhomogeneous cloud model in which 1% of the number density of acid drops are 
large (~13 ym) sulfur particles buried in a cloud just below that from which 
most of the . sulfuric acid scattering is observed, that is, at , T > 4 for 
visual wavelengths. Thus, these sulfur particles are well hidden in the visual 
and they would affect the spectrum primarily through their strong ultraviolet 
absorption. 

For the following reasons we suggest utilizing the infrared reflectivity 
of Venus to investigate such inhomogeneous models. First, scattering at 
intermediate backscatter angles by the cloud drops themselves is greatly 
reduced in the infrared absorption bands of H 2 SO 1 J. However, backscatter ing by 
nonabsorbing species, such as sulfur, is not so reduced and thus can be more 
easily detected. Second, the absorption is so strong that single scattering 
is the predominant radiative transfer process within the clouds. Third, the 
quite large size and irregular shape of the proposed sulfur particles will 
also assist their detection at 3 ym wavelengths. 

To evaluate our suggestion, we have calculated the reflectivity of Venus 
at 3.4 ym for two inhomogeneous models using a Mie scattering doubling routine 
which we modified to contain either an inhomogeneous drop with a sulfur core, or 
a bimodal particle-size distribution with large, nonspherical sulfur particles 
in the second peak. We have also calculated the homogeneous drop-homogeneous 
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cloud model of Hansen and Hovenier (ref. 3) for comparison. These models and 
representative calculations are presented in considerable detail in the next 
section in order to describe them and a new type of reflectivity function under 
a variety of conditions. In the third section we discuss the models; show the 
dependence of their backscattering on the core volume or particle area; evalu- 
ate the amount of sulfur that absolute observations at 3.4 ym could detect by 
comparison with them; and calculate the height to which the observations 
penetrate. 


MODELS AND DEMONSTRATION CALCULATIONS 


The Reflectivity Function 

We have chosen to define a reflectivity function, R(A,(f)), for the Venus 
cloud top that is an average of the directional reflectivity over the illumi- 
nated fraction of the planet's disk. This is a convenient choice because when 
absorption is strong, the shape of the phase angle dependence of the reflec- 
tivity function is very similar to that of the single scattering phase func- 
tion and so facilitates comparison of telescopic observations with cloud drop 
models. This definition has already been used by one of us (Pollack, ref. 4) 
to describe reflection spectra of Venus between 1.2 and 4.0 ym. 

The reflectivity of a cloud or atmosphere is fundamentally related to the 
wavelength of the incident light, X, and to the scattering angle, 6. Hence, 
the reflectivity function is an intrinsic function of x^avelength and angle. 

We write it as R(X,<j>) because the planetary phase angle <{> is the supplement' 
of 0 and is easily obtained from the Ephemeris. Unfortunately, the observed 
infrared flux from Venus depends both implicitly and explicitly on several 
other variables that mask the X, (j) dependence of the reflectivity. To 
elucidate that dependence, we write the monochromatic flux received at Earth, 
F(X, (j>, T, r^_^) , as the sum of two Lambertian surface brightnesses multiplied 
by the appropriate solid angles. The first term represents reflected light 
and the second term thermal emission 


f(X,<{>,T, r j ) 


Ra,<j.)F^_^(X,r^_^) 


+ cos j) 


irR^- 


IT 


+ 




c(X)W(X,T) 

ir 



In equation (1), R^ is the semidiameter of Venus, r^_e is the Venus-Earth 
distance, and is the Sun-Venus distance. The thermal spectral emit— 

tance is represented by a Planck function, e(X)W(X,T), in which a uniform 
effective temperature T is a good approximation for the entire disk of 
Venus on both the illuminated and dark sides. Equation (1) treats the spec- 
tral reflectance k(X, (j))Fg_^ (X,rg_„) in the same manner as the thermal spec- 
tral emittance after multiplication of the solar spectral irradiance, 

(XjTq^q) , by R(X,i|)) has converted the irradiance to a reflectance. Divid- 
ing both fluxes by tt sr converts them into Lambertian radiances for normal 
incidence. Thus, the reflectivity function R(X,ij)) is an average reflectivity 
over the illuminated portion of the Venus disk. 
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The phase-angle dependence of the atmosphere’s reflectivity is now con- 
tained solely in R(A ,<|))5 although also appears as (1 + cos in the 

solid-angle factor. That factor describes only the fraction of the area of 
the disk that is illuminated and has nothing to do with the optical properties 
of Venus' atmosphere. This point should be emphasized. Most phase functions 
of Venus are plots of reflected flux versus phase angle; therefore, they are 
heavily weighted by (1 + cos 4')/2. We have separated out that factor, so the 
reflectivity function just defined pertains only to an average segment of the 
illuminated atmosphere and is independent of illuminated area. 

For ease in comparing to related quantities defined by others, 

equation ( 1 ) may be rewritten as 




F(X,(f.,T,r5_e)S^2 ^ _ e(x)w(l,T) 


1 + cos 4 > 


9-e- 

^©—9 (X,r0_^) 


( 2 ) 


where \ " Ro/j^o-©* thermal term is insignificant at X < 3 ym. 

* * 

In this article we will be concerned with observations at X = 3.4 ym, where 
the reflected flux still dominates the thermal contribution. However, near 
inferior conjunction it may become necessary to correct observational data for 
'the thermal component when determining R(X,(fi). 

We now compare R(X,<j)) to the geometric albedo, p(X,<}>). In reference 5 
the albedo of Venus at unit distance from both the Sun and Earth is calcu- 
lated as 

log^Q p(X,(f),l) = 0.4[mg - m(l,(|))] - 2 log^gCsin oj) (3) 


where m^^^ is the apparent magnitude of the Sun at wavelength X; m(l,^) is 
the apparent magnitude of the planet at wavelength X if the planet were at a 
distance of 1 AU (astronomical unit) from both the Sun and the Earth ; and 
is the radius of the circle having the same projected area on the sky as the 
planet's disk at 1 AU from Earth. Converting magnitude to flux and uj to 
solid angle, the geometric albedo is 


p(X,<j),l) 


F(X,<}>,1 AU) 


892(1 AU)Fg,(X,l AU) 



(4) 


where Fg,(X,l AU) is the solar spectral irradiation at Earth and F(X,({),1 AU) 
is the flux that would be received from Venus if r^_® = 1 AU. But 


f©_9 (X , rg_^ ) 


= F^(X,1 AU) 



(5) 
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and 


F(X,(j),rg_e) ^ AU) 

S^^Cl AU) 


so that 


p(X,({)) 


F(X,<i>,r5_e) 

^9 (^^-© )^o— ^ > ^0— ^ ) 


( 6 ) 


(7) 


This definition is valid at wavelengths where thermal emission is negligible. 
Where that is the case we can omit the thermal term in equation (2) , and then 

- ^1 t coJ A 

Thus, our reflectivity function is just the geometric albedo divided by the 
illuminated fraction of the planetary disk. 

The Bond (or spherical) albedo of a planet is the probability that a 
photon will be scattered by the planet . It is defined as the ratio of the 
flux reflected at all angles to the incident solar flux at (|) = 0, and 
expressed in reference 5 as 


A(X) = p(X,0)q(A) 
where q(X) is the phase integral 


r 


q(X) = 2 j $(X,<())sin((f>)d({) 


(9) 


( 10 ) 


In equation (10), $(X,({)) is the flux density seen by the observer normalized 
to unity at full phase; reference 5 expresses this in magnitudes as 


loglO ^(X,^) = 0.4[m(X,l AU,0) - m(X,l AU,<J))] 
In flux units this is 


(11) 


$(X,<}.) 


F(X,j>,l AU) 
F(X,0,1 AU) 


which, from equations (4) and (8) is simply 


3>(X,ii.) 


pg.o *) 

p(x,0) R(x,0) 


( 12 ) 


(13) 
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Thus, equation (9) for the Bond albedo can be written in terms of the reflec- 
tivity function as 


A(X) = 2R(X,0) 




1 + cos 


&a,o) 


sin cj) d<i) 


•'0 


• or 


A(X) 


j R(X,<}))(1-+ cos ^)sin c|) dif* 

*'o 


(14) 


The Bond albedo of Venus is difficult to evaluate from observation alone 
because at very small phase angles the planet passes behind the Sun. Calcula- 
tion of R(X,<|>) from theoretical models can then be used as 4' approaches 
zero to supplement observations in the evaluation of A(X) through (14). 

In summary, we have defined a reflectivity function that is an -average 
directional reflectivity over the, illuminated part of the Venus disk and whose 
shape is similar to that of the single scattering phase function for strong 
absorption. . It is not the usually defined planetary phase function, and- its 
shape is very different from that of the geometric albedo. It can be easily 
integrated to produce the Bond (or spherical) albedo. 


' Homogeneous Drop-Homogeneous Cloud Model 

In recent years our understanding of the Venus clouds has been vastly 
improved through calculation of the multiple scattering of polarized light in 
a homogeneous model atmosphere (ref. 3).. To provide a framework in which to 
place the modifications that produce the Inhomogeneous models which follow, we 
first outline the homogeneous calculation as described by Hansen and Travis 
(ref. 6). This calculation is usually divided into two parts-; (1) single ■ 
scattering by independent scatterers within a small volume element, and 
(2) multiple scattering by homogeneous layers of the entire atmosphere. 

Because our Inhomogeneous modifications affect primarily the single scattering 
part of the calculation, we will describe that part quite fully. 

At the heart of the single scattering problem lie the Mie electric and 
magnetic multipole coefficients, ajj and bn-- They are formed from combinations 
;of Ricatti-Bessel functions, 'and are functions of only the complex index of 
refraction, m = n^. - in^, and the size parameter, x = 2iTr/X, where r is the 
particle radius and X is the wavelength. The scattering and absorption 
cross sections, as well as the single scattering phase function, can be found 
from these coefficients following the Mie theory (see refs. 7, 8). 

Let the single scattering cross section Og be written as 

( 15 ) 
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where Qg is the efficiency factor for scattering. Similar expressions can 
be written for the absorption and extinction (scattering + absorption) cross 
sections. Van de Hulst (ref. 7) has shown that 


Qs = 


X" 


CO 

J^(2n + DCja^P + |b^|2) 


(16) 


n=l 


and 


OO 

Qe = ^ S (2n + l)(He(an + b^) (17) 

^ n=l 

where Qg is the efficiency factor for extinction. The efficiency- factor for 
absorption Qg can be obtained from Qg - Qg, Van de Hulst also writes the 
single-scattering phase function P^^(r,6) for the unpolarized intensity as 

'’(r.9) ■ [|Sl<r,6)P+ ISjtr.e)^] (18) 

where 


- 2 n(n + i) f^n^n^''°s + VnC^os 0)] (19) 


S2(r,0) - 2 nS + 1)' + an'^n^^os 6)] (20) 

n=l 

The superscript 11 indicates that P^^(r,6) is the first element of the 4 x 4 
phase matrix which completely describes the angular distribution and polariza- 
tion of the scattered light. Here iTg and are functions of 0, the scat- 
tering angle, and can be written in terms of the Legendre polynomials. Note 
that 0 is the supplement of the planetary phase angle 4> used in the pre- 
ceding section. 

In practice, we used a computer program due to. Dave, ^ which computes the 
Legendre polynomials and Ricatti-Bessel functions from recurrence relations. 
The coefficients iTg and are calculated for specified scattering angles 
using straightforward upward recurrence relations based on the well-known 
properties of the Legendre polynomials and their derivatives. 


Dave, J, V. : Subroutines for Computing the Parameters of the Electro- 

magnetic Radiation Scattered by a Sphere. IBM Report 320-3237, 1968. 
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For specific particle radii, the Mie coefficients are calculated from: 


an(z) 


bn(z) 




n n— 1 

. m ^ xJ 

II 

1 

L m X. 



l^mAj^(z) + <ReC^(x) - (ReC^_^(x) 
[mA„(z) +l] e^(x) - ?^_^(x) 


( 21 ) 


( 22 ) 


where Ajj(z) = <J>^(z) /4>n^(z) , z = mx is the complex size parameter, and 
i|jj^(z), Ricatti-Bessel functions. Here A^, the logarithmic deriva- 
tive of Been used to avoid a computer overflow problem. Unfor- 

tunately, calculation of that derivative by means of an upward recurrence 
relation' has been found by- Kattawer and Plass (ref. 9) to involve a numerical 
instability (growing errors) for large values of n/z. They suggested a sepa- 
rate downward recurrence calculation of the derivative, which we use, to 
avoid the instability. 


The single scattering problem is completed by Integrating the phase func- 
tion and efficiency factors calculated for specific particle radii over the 
distribution of particle sizes expected in each volume element. Thus 


P 


11 

( 6 ) 


_AL 

2TTks 


r 


max 


|Sj(r,6) 


min 


S 2 (r, 0 ) [2]n(r)dr 


(23) 


and 


k 


s 

e 


^’^max 

’^min 


irr^Qg(r)n(r)dr 

e 


^’■max 

’^min 


Og (r)n(r)dr 
e 


(24) 


where n(r)dr is the number density of particles with radii between r and 
r + dr., and kg and kg are the size- integrated scattering and extinction 
coefficients. In the homogeneous model, we use the "standard" size distribu- 
tion employed by Hansen (ref. 10), 


n(r) = N 


(l-3b)/b -r/ab 
r e 

(ab)^^“^^^^\[(l - 2b) /b] 


(25) 


where N = total number of particles per unit volume, 1 = the gamma function, 
and 
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a = scattering effective radius ~ 


i 

1 


max 


rTTr^nCr)dr 


(25) 


max 


'rrr^n(r)dr 


'■mxn 


f 


max 


(r - a)^irr2n(r)dr 


b = scattering effective variance =■ 


mxn 


f 


(27) 


max 

I irr^n(r)dr 

“min 


Because the /n(r)dr may have unusual dimensions, depending on the form of 
n(r ) , it is useful to normalize the phase function as we have done in equa- 
tion (23) by utilizing the size- integrated scattering coefficient. Through 
this normalization a dimensionless phase function is retained regardless of 
the exact dimensions of J*n(r)dr. 

The multiple-scattering part of the calculation is based upon the doubling 
or adding method, the basis of which Hansen and Travis (ref. 6) have traced 
back as far as Stokes (ref. 11). The basic principle is simply that the 
reflection and transmission of a combination of two layers can be obtained by 
calculating successive reflections between layers if the reflection and trans- 
mission of each layer are known independently. The initial layers are taken 
to be plane-parallel and quite thin, with position inside them specified by an 

r" 

optical -depth, t kgn(h')dh’, which must be summed over all species in a 

unit volume. If the layers are of similar composition, the results for a 
thick homogeneous layer can be quickly built up in a geometric (doubling) 
manner. If some, or all, of the layers are dissimilar the results for a thick, 
inhomogeneous layer are obtained more slowly by adding the contributions from 
each of the dissimilar layers . A recipe for the doubling routine has been 
given by Hansen and Travis (ref. 6), and we use a similar version developed by 
one of us (Pollack, ref. 12). 

Calculation .of the light reflected by a planet ' s atmosphere is completed 
by integrating the multiple scattering from locally plane-parallel slabs of 
thickness x, such as those described above, over the illuminated portion of 
the spherical surface of the planet. The integration is carried out from the 
illuminated limb to the visible terminator by double Gaussian quadrature. 

The reflected flux thus calculated is divided by the product 

/l + cos (p \ / solar flux incident\ 

\ 2 /\on the planet / 


9 



to produce the absolute reflectivity function R(X,(()) defined in the previous 
section. 

The model calculation can be demonstrated in an interesting way by showing 
the strong similarity (mentioned earlier) between the shape of the single scat- 
tering phase function and that of the reflectivity function when there is 
strong absorption. In figure 1 we compare calculations of both at 0.55 ym for 
such a case (hypothetical Case No. 4 of Hansen and Travis, ref. 6). We also 
. contrast both of them in that figure with the geometric albedo, p(ij)). Here, 
and throughout the rest of this paper, we shall use the term backscattering to 
refer to light scattered at phase angles (f> < 70°, and the term forward scat- 
tering to refer to light scattered at phase angles > 110°. The similar 
shape of the two curves is quite reasonable for the situation where absorption 
limits the interaction of light with the clouds to only one or two scatterings. 
Both functions have a characteristically steep slope because absorption within 
the particle strongly attenuates backscattered light while fon^ard scattered 
light, due principally to diffraction, is not so attenuated. In spite of the 
absorption, a trace of a primary rainbow is indicated by the small peak at 
<}> = 18°. The different shape of the geometric albedo contrasts sharply with 
these curves indicating that R(X,(j)) is physically the more meaningful 
function. 

The definition of the single scattering phase function, equation (23) for 
P^^(0), indicates that only five independent parameters other than ({> or 0 
need to be specified- They are X, n^., n^, a, and b. The optical constants 
of sulfuric acid at various concentrations have been measured by Palmer and 
Williams (ref. 13) between 0.36 and 25 ym and by Jones (ref. 14) from 20 to 
50 ym. An 85% aqueous solution is currently indicated to be the correct con- 
centration of the Venus droplets (Pollack et al. , ref. 4). The effective 
radius and variance of the principal droplet-size distribution, a and b, have 
been well determined by visible polarization measurements (Hansen and Hovenier, 
ref. 3). As a demonstration of the success of the multiple scattering model, 
we compare it in figure 2 with the observations of Irvine et al. (ref. 5) at 
1.0635 ym where there is still a little H2S0jj absorption. The calculated 
reflectivity function has been multiplied by 0.899 to obtain the best fit to 
the observations. This factor indicates that a slight reduction of the 
single scattering albedo calculated from the optical constants of H 2 SO 1 J is 
required, most probably the result of our neglect of faint absorption by the 
very weak 1.065-ym CO 2 band (Boese et al., ref. 15). In figures 3 and 4 we 
demonstrate the high sensitivity of the model to variation of the real, and 
imaginary parts of the index of refraction in the H 2 S 0 [j absorption band at 
3.65 ym. 

Both polarized visible light and unpolarized 3-4 ym light reflected from 
the Venus cloud tops are particularly sensitive to the high level particle 
size distribution because they both result from primarily single scattering 
processes.’ (The degree of polarization is reduced by multiple scattering.) 
Hansen and Hovenier (ref. 3) have already shown that the visible polarization 
is very sensitive to the value of the effective radius a and is fairly sen- 
sitive to that of the variance b. The value they determined for b (0.07) is 
. small by terrestrial standards and could have interesting implications for the 
dynamics of the Venus atmosphere. 
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In figure 5 we demonstrate the sensitivity of the homogeneous drop- 
homogeneous cloud calculation at 3.65 ym to variation of the effective radius. 
The transition from strong forward scattering by large drops of 2.7 ym. effec- 
tive radius to the Rayleigh scatter by the 0.3 ym drops is quite apparent. 

The region of backseat tered light (J) < 70° shows high sensitivity to the 
value of .a. In figure 6 we demonstrate the model's dependence on the effec- 
tive variance. Unfortunately, it appears to be so small that Hansen and 
Hovenier' s (ref. 3) determination of b will not be well tested in the 3-4 ym 
absorption band. However, their value for the effective radius clearly will 
be. 


finally, we contrast the 3,4 ym reflectivity function with the visual one. 
In figure 7 we plot them both using a = 1.05 ym, b = 0.07, 
m(3.4 ym) = 1.338 - 0.1761, and m(0.55 ym) = 1.44 - 0.0000241. The vastly 
different shape and level of the two curves is striking. In particular, the 
hack scatter by the homogeneous model is more than 40 times smaller at 3.4 ym 
that at 0.55. Such a large reduction of the signal from the acid drops opens 
up the possibility of detecting small backscattering contributions by aerosol 
species which do not absorp in the infrared. 


Sulfur: Its Solubility and Optical Constants 

As the remainder of this paper hinges upon the concept that sulfur does 
not dissolve in sulfuric acid (personal communication, 0. B. Toon, Ames 
Research Center, 1978), we searched the International Critical Tables (ref. 16) 
and other sources for a quantitative measure of its solubility In the acid. 

None was found, but Fike (ref. 17) states that the sulfur coatings of 
concentrated-acid containers were generally unaffected after several years of 
service. We are indebted to Dr. T. W, Scattergood of SUNY-Stonybrook, who 
performed the following experiment and observations for us. 

An amount of chemically pure, powdered sulfur (sublimed) equal to 2.8 mg 
was dropped into 250 ml of aqueous solution which was 88 % H 2 SO 4 by weight. 

The sulfur, despite its slightly greater density, stayed on top of the solution 
until it was dispersed throughout by mechanical stirring. The sulfur did not 
appear to be readily wetted by the acid. In about 1 day much of it had settled 
to the bottom of the beaker although some still remained on the surface. After 
15 days at room temperature the distribution of sulfur powder was about the 
same, and less than about one fifth of it by visual estimate had disappeared. 
Thus an upper limit on. the mass of the sulfur that dissolved in the 88 % acid 
solution during 15 days was 1.3 parts 'per million. This limit may be uncer- 
tain by a factor of 2 because of the method of visual estimation employed. 

This experiment, together x>?ith the observation by Fike (ref. 17), provides a 
factual basis for the calculations of inhomogeneous models of the Venus clouds 
which follow. 

The optical constants of sulfur in the infrared are required to evaluate 
these inhomogeneous models, so we now describe how they were obtained. As 
will be shown, the dependence on their exact value is not strong. 
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Although the constants are well determined at 0.5893 ym (ref. 18, 
p. E-219) they are relatively unknown at other wavelengths. However, spectra 
of orthorhombic sulfur, Sq, have been obtained from 0.1 to 200 ym by several 
authors. In particular, spectra from Nyquist and Kagel (ref. 19), Khare and 
Sagan (ref. 20), the Thermophysical Properties of Matter (ref. 21), and 
MacNeil (ref. 22) are free of all absorption features (other than those of 
expected contaminants) between 2 and 7 ym, and the nearest strong sulfur 
absorption is at 11.8 ym. Thus, we feel confident that the imaginary index 
near 3.4 ym is quite small. MacNeil estimated from her transmission spectra 
that the exponential absorption coefficient a, between 2 and 10 ym, was about 
4 cm~^. Hence, we use 



with X in ym. 


For the real part of the refractive index, normal dispersion can be 
assumed because absorption has been shown to be negligible. As we are also- 
far removed from strong absorption features, we need employ only the first two 
terms of Sellmeyer’s equation (ref. 23, p. 96) to calculate the real index to 
the accuracy we require. Thus we use 


n_^ = 1 + 


BX^ 

X2 - Xq2 


(28) 


with B determined from the known index at 0.5893 ym, and the fundamental 
wavelength Xq = 0.13 ym from Cook and Spear (ref. 24). This gives n^ = 1.96 
at 3.4 ym. MacNeil (ref. 22) calculated a reflectivity of 0.075 between 2 and 
10 ym for one surface of a single crystal, which requires n^ = 1.75 to 
satisfy the equation for reflection at normal incidence [ (n - l)2/(n + 1)2], 
Noting that the dielectric constant of orthorhombic sulfur = 4.0 (ref. 18), 
we select a value for n.^. of 1.80 at 3.4 ym. Young (ref. 1) states that the 
real index of his sulfur particles is 1.9 (presumably at 0.9 ym, the wave- 
length of the Venera 9 and 10 nephelometers) . 


Inhomogeneous-Drop Model 

The Mie theory for single scattering by an inhomogeneous drop was first 
worked out by Aden and Kerker (ref. 25) and by Guttler (ref. 26). Recently, 
Wickramaslnghe (ref. 27) has reformulated the equations for inhomogeneous Mie 
coefficients in a more explicit manner, although we note that a printing error 
has nearly obliterated a minus sign before the first term of the second line 
of his equation (3.45). As Kerker (ref. 8) noted, the variety of numerical 
calculations that might be explored with the inhomogeneous drop model is 
immense because the size parameter and refractive index of both the core and 
coating are variable parameters. Although he considered many examples, Kerker 
did not treat the one applicable to the suggestion of sulfur in the Venus 
clouds . 
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Consider an inhomogeneous drop (fig. 8 ) containing a spherical, homoge- 
neous core of radius Rg and complex index m^ j surrounded by a concentric 
coating of outer radius R and index m 2 » For the example of sulfur coated 
by H 2 SO 4 this model describes at 3. 4 ym a transparent dielectric core sur- 
rounded by a strongly absorbing dielectric coating. In the specific case of 
Venus, the visible polarization observations set the effective outer radius R 
at about 1 ym. -Thus, only the core radius remains as a variable 

parameter. 

There are at least four different ways, however, in which Rg can vary 
relative to R as the integration over the particle size distribution is 
performed. In most cases, the core-size dependence on R can be written as 

Rg = -1 < p (29) 

where C is a constant. In the first case, p = -1 so that R^ remains 
constant throughout the size distribution; particles in the distribution 
smaller than C are entirely homogeneous core material, with coating applied 
only to drops larger than C. In the second case, p = 0, so that '• R^ scales 
as a fixed fraction of R throughout the distribution. We call this the 
scaled-core model. In a third case, p > 0, hence Rg grows faster than R 
across the size distribution and will eventually equal R, say at some value 
V. In this case, C = V f*, and particles larger than V are entirely core 
material; we call this the growing-core model. A fourth case, which we shall 
call the constant coating model, has a constant thickness t of liquid 
coating the core in those drops large enough to contain a core, that is, 
where R > t. For R < t, the drops are entirely liquid material. To describe 
this model, equation (29) is modified slightly so that Rg = R - t. We note 
that for small values of t this model has the same core-size dependence on 
R as does the scaled-core model, that is p = 0 and C = 1. 

New Mie coefficients must be employed to describe the inhomogeneous drop. 
Once they have been obtained, the calculation proceeds as in the homogeneous 
model. Equations for the required coefficients are given, for instance, by 
Kerker (ref. 8 , pp, 189-195) and Wickramasinghe (ref. 27, pp. 29-31). Rather 
than reproduce them here, suffice it to say that we have evaluated them for 
the first three orders. 

We have not attempted to derive a general recurrence relation for the 
inhomogeneous Mie coefficients because in most cases they need to be’ used only 
in the first few terms of the series for the phase function and efficiency 
factors. The remaining terms can be obtained in the usual manner from 
homogeneous coefficients. The reason for this follows from the concept of ray 
localization for x > 1. According to this principle, the nth term in each 
series corresponds to the contribution from the light ray that passes the drop 
at 'a distance, mr, of nA/2ir from its center (fig. 8 ). Thus, terms where n 
is less than 2 Trm 2 Rg/A are strongly influenced by the core material and must 
be calculated from inhomogeneous Mie coefficients. ■ Higher order terms corres- 
pond to rays that miss the core; they can be calculated from the homogeneous 
model. Thus it will be useful to define an inhomogeneous drop parameter, nl, 
to indicate how many Inhomogeneous coefficients are required. For the specific 
case of Venus at 3.4 ym. 
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when using Hansen and Hovenier's effective drop radius of 1.05 pm (ref. 3). 

We shall shortly show that only the first three terms need to be calculated 
from inhomogeneous Mie coefficients for the narrow particle size distribution 
of Venus. 

In the more general case, the value of nl will depend upon the manner in 
which Rq scales with R, and also on the width of the particle size distri- 
bution. However, by defining nl we avoid the difficulty calculating 

inhomogeneous coefficients described by Kerker (ref. 8, p. 200). In that 
situation, n was set by the value of R rather than by R^, which caused 
n >>® ZQore those Riccati-Bessel functions and derivatives associated with 
the core as Rq 0. As mentioned earlier, Kattawar and Plass (ref. 9) found 
that calculations employing upward recurrence relations are numerically 
unstable in such a case. Thus, by using inhomogeneous coefficients obtained 
through upward recurrence only so far as nl, and then using homogeneous 
coefficients obtained from downward recurrence for the remaining terms, the 
numerical instability is avoided for most values of Rq/R. A thin film coat- 
ing on particles with a broad size distribution is an extreme example where 
it would be desirable to have a general recurrence relation for the inhomoge- 
neous coefficients because of the large value of nl. 

The multiple scattering Mie program for planetary reflection (described 
earlier) was modified to calculate inhomogeneous Mie coefficients through the 
third order following Wickramasinghe (ref. 27). It was also, modified to 
permit Rq to vary according to equation (29). The inhomogeneous drop pro- 
gram thus created was tested by requiring it to produce the same results as 
the homogeneous drop calculation to five or more significant figures in four 
equivalent test cases. The first equivalence test had mj^ = m2 = the index of 
H2S0i^, and R^ had a constant value of one half the modal drop radius (or 
0.415 pm). The second test case had m^ = the index of H2SOJJ, m2 = 1, and Rq 
had a constant value of 4 pm (the size distribution peaks at 0.8 pm). The 
third test case had mj = sulfur, m2 = H2SO4, and a scaled core where 
Rq = 0.01 R. It was for such very small cores that Kerker (ref. 8, p., 200) 
experienced difficulty with upward recurrence as mentioned above. We also 
calculated his specific case (X = 2.5 pm, m^ = 2.105-0, m2 = 1.482-0, x = 10, 
Rq/R = 0.005, b = 0.07) and found no difference from our corresponding homoge- 
neous calculation. In the fourth test, m^ = sulfur, m2 = H2SO4, and_the con- 
stant coating thickness was 2 pm, a value so large that less than 10~^ of the 
drops had cores. 

In a fifth equivalence test the scaled core size was as large as the drop: 
Rq == R, mj = H2SO4, and m2 = 1. This case thoroughly tests the inhomogeneous 
drop program because three different values of z remain defined throughout 
the calculation and because |z(,Qj.g| approaches 5 at the large end of the size 
distribution, thus fully .utilizing the higher-order inhomogeneous Mie coeffi- 
cients. The result of using more inhomogeneous coefficients in this test is 
clearly shown by the change of the single scattering phase function in fig- 
ure 9. Other functions, such as the anisotropy parameter, the single-scattering 
cross section, and the albedo similarly approach the homogeneous result as more 
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inhomogeneous coefficients are used. Inasmuch as this is an extreme test 
(Rq = R) and the residual deviation is less than 2% after the third inhomoge- 
neous term, we consider nl = 3 more than adequate for the Venus calculation. 

Finally the inhomogeneous drop program was tested by using it to recal- 
culate the scattering and extinction efficiency factors of inhomogeneous drops 
published by two other authors. Since theirs were monodisperse calculations, 
our program was temporarily modified to bypass the particle size integration. 

In each test when the core radius became so large that nl > 3, the test com- 
parison failed significantly, as it should. Espenscheid et al. (ref. 28) 
calculated a nonabsorbing case where m2 = 1.481, m^ = 2.105, and x = 10. 

They used upward recursion wherever possible (small values of n/z) which may 
account for the small, but consistent difference of 2/3% between their values 
of Qscat Rq/®- approached zero. Fenn and Oser (ref. 29) calcu- 

lated the inverse case to ours, that of an absorbing core (m^ = 1.59-0.661) 
covered by a nonabsorbing coat (m2 = 1.33). At values of Rq/R equal to 0.5 
and 0.833, with x between 1 and 4, our calculations of Qg and Qg agreed 
with theirs to within our ability to read their graphs (about 4% or better). 

They also calculated absorption efficiency factors for an absorbing coating 
(m2 = 1.473-0.1691) with the same absorbing core. At x = 12.3 our values 
of (obtained from the difference between our calculated values of and 

Qg) at Rq/R = 0.0273 and 0.0741 agreed with theirs to 5 significant figures. 
However, comparison with certain plotted data in Wickramasinghe (ref. 27) and 
Pilat (ref. 30) produced agreement only on the order of 10%. 

We now want to demonstrate the dependence of the reflectivity function on 
the variable parameters of the inhomogeneous drop. I-ihere the homogeneous drop 
has five independent parameters, the inhomogeneous drop has eight, as well as 
the manner in which the core radius varies over the size distribution. For 
cores that do not occupy the entire drop, it is reasonable to assume that the 
reflectivity function of the inhomogeneous drop model varies with wavelength, 
effective radius, variance, and coating index in about the same way as it does 
for the homogeneous drop; thus we will not demonstrate the dependence on those 
parameters. Variation of the three additional parameters of the inhomogeneous 
drop — Rq, and the real and imaginary parts of the core index — can produce 
a dramatic effect on. the reflectivity function depending on their value. How- 
ever, for cores of sulfur the imaginary index is so small (^10“*^) that varia- 
tion of this parameter over a reasonable range has no effect whatever. The 
real index of sulfur is larger than that of the H2SO4 coating, so its varia- 
tion could be significant. In figure 10 the real part of the core index is 
varied over the range of possible values (1.75-1.95), It can be seen there 
that this range of variation does have a noticeable, but small, effect on the 
reflectivity. The model's sensitivity to this core parameter is much less 
than it is to the real index of the coating (fig. 3) . 

We now compare the four previously described ways of varying the core 
radius during integration over the size distribution. This is done in figure 11 
for moderately large cores. Principally because of their smaller volume of 
absorbing material, the scaled- core and constant-coating models have signifi- 
cantly greater reflectivity than either of the other tx-70 models at most phase 
angles. We noted earlier that for small coating thicknesses both of these 
models have the same- core-size dependence on drop radius (p = 0); thus, it is 
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■not surprising that their reflectivities have a similar dependence on phase 
angle. It is a bit surprising, however, that the magnitude of the reflectivity- 
calculated for these two models is so similar (see also figs. 12 and 13). The 
constant in equation (29) is a measure of core size for the scaled core model 
(= Rq/R). In the constant coating model Rq/R = 1 - t/R = 1 ~ t/a, where we 
have used R = a because of the small width of the size distribution. By 
comparing figures 12 and 13 it can be seen that as the core increases in size, 
the two reflectivity functions become nearly equal, and furthermore the two 
measures of core size also approach equality at C = 1. Thus, for large cores, 
not only do the txi/o models have the same core-size dependence on drop radius, 
they effectively have the same core size. The narrowness of the size distri- 
bution helps produce this effect. The resulting near-equality of the reflec- 
tivity functions of the two models is thus not so surprising. 

The constant-coating, the scaled-core, and the constant-core models 
should also be compared (figs. 12, 13, and 14) at the opposite range core- 
size variation as Rq/R approaches zero. Both the constant-core and scaled- 
core models regularly approach the homogeneous-model reflectivity function as 
core size diminishes. However, the reflectivity of the constant-coating model 
at small phase angles dips somewhat below that of the homogeneous and scaled- 
core models when Rq/R is <0.33. This small enhanced absorption occurs when 
thei constant coating thickness is approximately one-fourth the wavelength. 

The growing core in figure 11 is described by R^ = R^/(2 ym) , so that 
Rq R at 2 ym. For convenience in doing this particular calculation, the 
sj-ze. distribution was cut off at 2 ym where its value is more than three 
qrders of magnitude below the peak value. We note that the reflectivity of 
the growing core model has a steeper slope than the other three models because 
the larger drops become relatively less absorbing than the smaller ones. The 
reflectivity of the constant core example has an intermediate slope. Although 
these differences are noticeable, they are much smaller than those due to 
changing core size. 

Regardless of the mode of variation of the core radius, the reflectivity 
function is most sensitive simply to the effective size of the core itself. 

This strong sensitivity is clearly shown in figures 12-14, and is based on the 
greater volume of highly reflecting, nonabsorbing material in the larger cores. 
When Rq/R is greater than 0.5, this sensitivity is particularly strong, and 
the reflectivities of the scaled-core and constant— coating models are nearly 
equal, as noted above. In those models, both the level and the slope of the 
reflectivity function change rapidly for R^/R >0.7. Finally, drops with 
cores < 0.1 R are nearly indistinguishable from homogeneous particles, with 
the small variation of the constant-coating model noted before being, the only 
exception. 


Inhomogeneous-Cloud Model 

If large particles of elemental sulfur are hidden at large visual optical 
depths as suggested by Young (ref. 1), then they might be detected in the 
i^f^sted because a small percentage of them could be brought up to the level 
of the 3.4 ym observations by convective overshoot or other turbulent 
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processes. There the particles would be particularly noticeable because 
absorption in the acid drops has already reduced the amount of scattered 
light. A small contribution to the scattering by nonabsorbing particles should 
be easily detected, particularly so in the intermediate backscatter directions 
where the droplet attenuation is strongest. Because both the large size 
(~13 pm) and the irregular shape of the postulated sulfur particles enhance 
the infrared backscatter relative to that of acid drops, these characteristics 
would further assist sulfur particle detection. 

Since most infrared observers do not resolve planetary detail, we must 
consider the postulated sulfur particles to have a uniform horizontal distri- 
bution in spite of the known ultraviolet contrasts. We also assume that they 
have a uniform vertical distribution within the top few optical depths for 
reasons of mathematical convenience. Hence, in this first approximation we 
consider a uniform cloud of proposed sulfur particles that is well mixed with 
the H 2 SO 4 cloud drops. A bimodal particle size distribution describes this 
situation quite well because the inhomogeneity in particle size, indices of 
refraction, and number density have at least as large an effect upon the 
planet's backscatter as does its spatial inhomogeneity. Thus we modify the 
homogeneous model described earlier to include sulfur simply by replacing the 
existing unimodal size distribution with a bimodal one that has H 2 SO 4 drops in 
the main peak near 1 pm and irregularly shaped sulfur particj.es in a second, 
smaller peak near 12- pm. A more general treatment in i^hich the size distribu- 
tions overlap is also possible, but is not called for by either the quality 
of current data or the hypothesis of Young's model. 

To modify the homogeneous model to incorporate sulfur particles we 
replace equation (25) by 

n^(r) = F^n(r,aj,b^) + F 2 n(r,a 2 >b 2 ) (31) 

where n(r,a,b) is given by equation (25) and Fj + F 2 = 1. Here F^ is the 
fraction of the total number of particles (acid drops plus sulfur) in the 
sulfuric acid peak at a^, and F 2 is the fraction of all particles, in the 
sulfur peak at a 2 - The effective radii and effective variances aj, 32 , bj, 
and b 2 are defined by equations (26) and (27) for the two different modes 
by considering the limits of integration, and to separate the par- 

ticles into two complete and exclusive distributions. In practice, however, 
a^, 32 , bj, and b 2 are variable parameters. 

Following equations (23) and (24), the size— Integrated phase function 
and cross sections become 
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(The numbers 1 and 2 within the brackets of equation (32) do not refer to the 
bimodal distribution; see eqs. (19) and (20).) Here Tq is an intermediate 
integration limit where we change both the particle shape and the index of 
refraction to values appropriate to those of solid, irregular sulfur particles. 
As already mentioned, our treatment of the bimodal distribution as completely 
separable is a simplification allowed by the state of current theory and data. 
It will be shown later that this scattering solution is quite insensitive to 
values of Tq between 3 and 7 pm. A more general treatment would merely 
require a different integral. 

At the temperature of the Venus cloud top (~240 K) , elemental sulfur can 
be expected to be frozen and thus have an irregular shape. The scattering 
behavior of irregular (nonspherical) particles of size somewhat larger than 
the wavelength differs systematically from that predicted -by Mie theory. In 
particular, a much smaller peak is seen in the directly backward direction 
(<j> = 0) and the scattering over intermediate backscatter angles is more uni- 
formly distributed. Cuzzi and Pollack (ref. 31) describe a semi-empirical 
method for obtaining tbe single scattering phase function and albedo of 
irregular particles given the Mie values for a set of equal volume spheres of 
the same refractive index. Their method uses a combination of physical and 
geometrical optics, and relies on comparison of parameterized phase functions 
with laboratory observations of the scattering by irregular particles of var- 
ious sizes and shapes. The semi-empirical theory fits these observations 
quite well with variation of only two parameters. The first is the size 
parameter, x^ = Zirr^/A, where the transition from Mie scattering by spheres to 
irregular particle scattering occurs. The second, FTB, is the slope of a 
straight line representing the angular dependence of that fraction of the 
scattered light that has been transmitted through the particle. Cuzzi and 
Pollack (ref. 31) have shown that values of these parameters characterizing a 
wide range of particle shapes do not vary so widely that "typical" values may 
not' he chosen. Thus, we have used their semi- empirical theory with 
Xq'= 7 (rQ =3.8 pm), a value characteristic of large particles, in our 
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calculation of the scattering by the sulfur mode of the bimodal distribution. 

In demonstrating and testing our calculations we have used FTB = 1-5, charac- 
teristic of cubes. However, a value of 3 is physically .more reasonable for 
large polyhedra, and we use that value later to set observational limits on 
other parameters. 

The modifications just described were tested by several reductions to 
known cases. First, using only a unimodal size distribution (a = 1.05 ym, 
b = 0.07), the spherical-drop part of the program was shorn to reduce to the 
H 2 SO 4 homogeneous drop, homogeneous cloud model reflectivity function at 1 ym 
(Case No. 3 of Hansen and Travis, ref. 6 ) and at 3.65 ym by setting very 

large (equal to 99). Still- using a unimodal size distribution, the irregular- 
particle part of the program was checked by requiring it to reproduce a pre- 
viously calculated phase function at A = 3.71 cm when the value of fell 

near the center of the range of integration. This calculation was related to 
laboratory observations of the scattering of microwaves by particles with size 
parameters similar to ours (fig. 5 of Cuzzi and Pollack, ref. 31). The 
bimodal distribution was tested in two separate cases by requiring it to 
reduce to the 3,65 pm absolute reflectivity function of the homogeneous model 
calculated earlier. In the first case, the second peak was simply eliminated 
(F 2 =0). In the second case, the first peak was eliminated (Fj = 0) and the 
second peak converted to H 2 SO 4 spheres (Xq = 99) located at 1.05 ym 
(a 2 = 1.05 ym) . In tests where previous calculations were reproduced, the 
agreement was good to five or more significant figures. In reproducing 
Case No. 3 of Hansen and Travis (ref. 6 ), the agreement was as good as our 
ability to read their figure 27 (about 3%). 

We will now discuss the dependence of the reflectivity function on the 
variable parameters of this inhomogeneous cloud model and point out those that 
are the most significant. Since the fraction of all particles in the second 
(sulfur) mode that we may expect to see at 3.4 ym will probably be less than 
the 1% suggested by Young (ref. 1), it is reasonable to assume that this model 
will have a dependence on the parameters of the sulfuric acid droplet mode 
similar to that of the homogeneous model. In addition to those five param- 
eters, the inhomogeneous cloud model presents seven more to consider. As with 
the inhomogeneous drop model, the dependence on the imaginary index of sulfur 
is negligible because its value is so small in the first place. The depen- 
"dence of the inhomogeneous cloud model on the real index of sulfur (fig. 15) 
was also found to be negligible. It was mentioned earlier that because the 
modes of the size distribution were separate, the dependence on Xq was small 
for Tq in the range between 3 and 7 ym. Of course outside that range, the 
sensitivity to Xq can be appreciable because Inappropriate shape and index 
changes are being made on the shoulder of one mode or the other. These effects 
can be seen in figure 16 where is varied from 1,6 to 11.5 ym. Unless 

stated otherwise, all succeeding calculations were made xvrith r^ = 3.8 ym 
(Xq = 7). The inhomogeneous cloud model is not as sensitive to FTB, the slope 
of the parameterized phase function of the transmitted light, as might have 
been expected. Figure 17 shows, however, that the dependence does occur where 
expected, at the largest backscatter angles, because the transmitted component 
dominates at large scattering angles when FTB is small. The models' dependence 
on the effective variance of the sulfur peak, shown in figure 18, is also 
small, Backscattering increases slightly with decreasing variance because the 
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"standard" size distribution employed here is asymmetric and narrowing the 
second mode creates proportionally more particles with larger radii. 

The backscatter i 70“) of the inhomogeneous cloud model is most sen- 
sitive to two parameters: the effective radius of the sulfur particles, a 2 , 

and the fraction of all particles in the sulfur mode, I'2- will be shown in 

the next section this strong sensitivity, depicted in figures 19 and 20, 
results from a nearly linear dependence of the model' s backscattered reflec- 
tivity on the total area of large particles. Thus curves with about the same 
reflectivity at <j) = 20“ are nearly identical, regardless of whether the 
varied parameter was &2 ^2* range of F 2 shown in figure 20 is from 

zero to 1/100, which spans the range from a homogeneous model to an inhomo- 
geneous model containing as much sulfur as was proposed by Young (ref. 1). It 
is clear from this figure that the backscattering of the inhomogeneous cloud 
model is sensitive to only a few percent of the proposed amount of sulfur. 


COMPARISON AND DISCUSSION OF MODELS 


Backscatter Dependence on the Sulfur Volume or Area 

In the preceding section we found that although there can be many param- 
eters associated with the different inhomogeneous models, there are only a- few 
that really govern the shape and level of the reflectivity function. We shall 
call these the principal parameters of each model; they are listed below. The 
first three' in each model refer to sulfuric acid drops. 


Homogeneous model 

a, n^, 

Inhomogeneous drop 


Inhomogeneous cloud 



The reflectivity function of the homogeneous model has a broad minimum 
near a phase angle of 30° (fig. 7); hence that is a region where scattering 
associated with the additional principal parameters of the inhomogeneous models 
will be most apparent. Accordingly, we investigate the dependence of the 
backscattered light on R^, a 2 j and F 2 at two backscattering phase angles: 

(J»g - 20° and ([ig = 40°. 

In figure 21 we have plotted the reflectivity at (|)g of both the 
scaled-core and constant-core inhomogeneous drop models versus core volume. 

We used a core volume 4/3ir (aRQ/R) ^ employing the effective radius a of the 
size distribution. Using the mean or modal radii increases the slope of the 
curves but does not otherwise change the results. It is clear from the figure 
that the backscattering of the constant-core model has a highly linear depen- 
dence on the sulfur c.ore volume which can be expressed as 

R(3.4, <j)g = 20°, 40°) = 0.0436 (Ro/r)3 + 0.00774 (34) 
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The difference between the backscatter at 20° and 40° is negligible so aver- 
ages have been used in equation (34). It can also be shovm that the back- 
scatter has a similar, but inverse dependence, on the volume of H2S0^,. Thus, 
we conclude that it is primarily the volume of absorbing coating that depresses 
the backscatter and causes the reflectivity function to be as steep as it is. 

It is quite apparent in figure 21 that backscattering by the scaled-core 
model is not linearly dependent on volume, having a much stronger dependence 
than that for > 0.4R. The const ant- coating model produces a, nearly identi- 
cal plot of R versus volume because for large cores it has the same depen- 
dence of core size on drop radius as does the scaled-core model, as noted 
earlier. The difference between the constant-core model and these other two 
models probably arises from the fact that the constant-core model has a con- 
stant volume of sulfur at each fractional size level, whereas the other two 
models have a variable volume of sulfur through the integration over particle 
size. Since the backscatter is proportional to sulfur volume, the size inte- 
gration permits the larger fractional cores to dominate the reflectivity of 
the scaled-core and constant-coating models. 

We noted earlier that the sensitivity of the inhomogeneous-cloud reflec- 
tivity function to both F 2 and a^ was due to its dependence on the total 
area of sulfur particles. The inhomogeneous cloud reflectivity function versus 
both parameters is plotted in figure 22 (the upper four curves). The nearly 
linear dependence on F 2 can be expressed by 

R(3.4, = 20° to 40°) = 20.8 F 2 + 0.00774 , with a 2 = 12 ym (35) 

and the dependence on a 2 (in microns) by 

R(3.4, (j)g = 20° to 40°) = 4.42xl0”5(a2)2 + 0.00774 , with F 2 = 3.2x10“** 

(36) 

Again, the difference in backscatter between 20° and 40° is so small that aver- 
ages have been used to obtain these equations. To determine these relations, 
we found that in the computer calculation the integration over particle size 
had to be carried to very large values (x = 85) to get all the backscattered 
light. As in previous calculations of this model, we have used FTB = 1.5 
for demonstration purposes. Because this may be physically unreal (large par- 
ticles do not often have square corners) , we have also calculated R as a 
function of F 2 with FTB = 4. These are the lowest two curves in figure 22. 
Increasing FTB reduces the slope of the dependence on F 2 , but does not 
remove the dependence or its linearity. These results are all consistent with 
the observation that non-absorbing, high-index particles that are considerably 
larger than the wavelength are good reflectors, and the reflection is propor- 
tional to the total area of the reflector. 


Comparison of Models 

Similarities- All three models are similar in that at 3.4 pm they have no 
glory or rainbows as <}> approaches zero. The homogeneous model produces only 
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a small monotonic increase of in R as <}> decreases from 40® to 0,5°. 

The inhomogeneous models have an even smaller increase. This result is con- 
sistent with most of the particles being considerably smaller than the wave- 
length. Both inhomogeneous models are similar in that variation of any of the 
"inhomogeneous" parameters changes the backseat tering more than the forward 
scattering because the forward scatter is due principally to diffraction by 
the abundant H 2 SO 1 J drops. Variation of in the inhomogeneous drop (up to 

values of about 0.6R) produces changes in the scattering that are somewhat 
similar to those produced by changing either n^ or of the acid in the 
homogeneous model. This happens because the sulfur volume determines the 
volume of acid remaining in the drop. Thus, small sulfur cores could have an 
effect on 3.4-ym observations similar to a reduced concentration of the acid. 
Fortunately the wavelength dependence of an H20-H2S0x^ solution is sufficiently 
different from that of an S-H 2 SO 11 mixture that large errors in the deduced 
concentration are not possible when observational spectra covering both liquid 
H 2 O and H 2 S 0 i^ bands are analyzed. 

Differences- Variation of the radius of the acid drops of the homogeneous 
model affects the 3.4 pm reflectivity function in a manner which is distinctly 
different from that due to varying other parameters of the same or different 
models. (Compare fig. 5 with figs. 14 and 20, for example.) Such changes are 
predictable. First, as the particles become small compared to A the R 
tends toward an isotropic distribution similar to Rayleigh scattering. Second, 
as the particles approach A in size, diffraction enhances forward scattering 
while absorption within the particles reduces the backscatter, resulting in a 
very steep reflectivity function. 

In figure 23 we contrast the effects of varying the other principal 
parameters of the three models. Variation of F 2 and a 2 of the inhomogeneous 
cloud model changes the shape of the reflectivity function through a purely 
backscattered enhancement; varying Rq of the sulfur core, or or of 
the acid, primarily affect the overall level of the curve. In summary the 
large, non-absorbing, high-index particles in the second mode of the inhomo- 
geneous cloud act roughly like an isotropic reflection added to the droplet 
scattering, while the volume of non-absorbing core inside the inhomogeneous 
drops affects the overall level of the reflectivity function by excluding the 
sulfuric acid absorber. 


The Detectability of Sulfur 

Minimum detectable amounts of sulfur can be set by comparison of inhomo- 
geneous reflectivity functions with that of the homogeneous acid drop. But to 
do that, the ratio of the mass density of elemental sulfur to that of a speci- 
fied concentration of sulfuric acid first must be evaluated. We will use an 
85% concentrated solution with a density of 1.71 gm/cm^, and note that the 
density falls by less than 5% if a 75% concentrated solution were used. The 
major uncertainty in the mass density ratio arises from the unknown molecular 
form of elemental sulfur on Venus. However, Meyer (ref. 32), lists in his 
table 3.1 the density of nine different allotropes of sulfur (ranging from 
^6 'to ^oo) and shows that they vary from 1.94 to 2.21 gm/cra^, a range of only 
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13%. Thus we pick a middle value of 2.1 gm/cm^ which produces a value for the 
S/H2S0£^ mass density ratio of 1.23 with an uncertainty of ±8%. 

Detectability of sulfur in the Venus clouds, in the sense that we shall 
use it, refers to. the ability to observe an enhancement of the 3.4 ym back- 
scattering relative to that calculated for a homogeneous drop. For example, 
figure 24 shows a band 25% higher than the homogeneous drop reflectivity 
function. Just within that band lie back scattered reflectivities of inhomo- 
geneous drop, models with scaled core radii as large as 0.38 R, and 
inhomogeneous cloud reflectivities with F2 = 0. 75><10~^. Thus, 3.4 ym obser- 
vations at ij) < 100“ with total (absolute plus internal) errors <25% could 
distinguish the homogeneous model from those inhomogeneous drop models whose 
scaled core radii were >0.38 R. Because 25% error bars placed upon those 
curves with Rq < 0.38 R would include the homogeneous curve, such observa- 
tions could not distinguish, hence could not detect, inhomogeneous drops with 
cores smaller than 0.38 R\ Similarly, figure 24 indicates that ±25% total- 
error observations at = 20“ to 40“ could distinguish the homogeneous 
model from inhomogeneous cloud reflectivities where F2 > 0.75x10”** (with 
a2 = 12 ym, FTB = 1.5). 


The ratio of the detectable mass of sulfur mg to that of H2SO4 solution 
“H2SO4 directly related to the minimum distinguishable values of F2 

or Rq/R set by the total error of the observations. This is particularly 
simple for the scaled core inhomogeneous drop model where the ratio of the 
total volume of sulfur (V^) in the atmosphere to the total H2S0t^ volume 
(Vh 2S0 £j^ is the same as it is in each drou because there is a core of fixed 
fractional size in every drop. Thus, 


V 


H2S0^ 



1 

(R/Rq)3 - 1 


(37) 


and so the minimum detectable mass ratio for inhomogeneous drops is 

1.23 


m. 


mid 




Ps^s 


drops 


Ph2S0Jh2S0j^ 


(R/R„) - 1 


(38) 


where we have used a density ratio PS/PH2SO4 ” 1-23 as discussed earlier. 
For the inhomogeneous cloud model, the minimum detectable mass ratio is 


m. 


A . E - 
mic m. 




Ps^l^s 
~ 3 

clouds ’^H2S0 ij ^1 ^H2S04 



( 39 ) 


since Ng/%2S04 " ^2 we consider to be specified by the work of 

Hansen and Hovenier (ref. 3).. 
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To take the case from figure 24 as an example, inserting R/Rq “ (0.38)“^ 
into equation (38) gives \a.xd 7.1% for the inhomogeneous drop. 

That is, 25% or better observations could tell inhomogeneous drops with more 
than 7,1% sulfur from homogeneous acid drops. More accurate evaluation of 
Rq could be made by referring to plots of R versus volume, such as 
figure 21, but that is not necessary at this point. 

The inhomogeneous cloud model has two unspecified parameters,' a .2 and R 2 * 
Young (ref, 1) has proposed that ^2 > lOa^^, so we will evaluate equation (39) 
at two values of ^ 2 "^ mentioned earlier, values of 

FTB = 1.5 may be unreal for large particles, so we use now a value of 3 
appropriate for large polyhedra. This is consistent with the shape of the 
double pyramid habit of orthorhombic sulfur shown by Adams and Spear, (ref. 33). 
Because of the high sensitivity of the back scattered reflectivity to F 2 , it 
is both impractical and inaccurate to evaluate the minimum distinguishable 
(detectable) value of F 2 from families of curves such as those in figure, 20. 
Since R is directly proportional to F 2 , the minimum distinguishable value 
can be extrapolated from a linear plot as shown in figure 25 for a 2 = 18 pm. 
Again we note that the difference between <p = 20° and 40° is small enough to 
permit averaging, so that if a 2 = 18 pm we find from the figure that' the min- 
imum distinguishable value of F 2 is 6.1x10"^ at the ±25% total-error level. 

It is of course proportionately smaller at the ±15% error level. 'Similar cal- 
culations for a 2 = 9 pm give the minimum distinguishable value of F 2 -as 
2,2x10”^ at the ±25% error level. Thus, depending on the size of the postu- 
lated sulfur particles we expect that their minimum number density detectable 
by ±25% total-error observations at is between 6x10“^ and 22x107^ times 

the number density of acid drops. This is about two orders of magnitude 
smaller than that proposed by Young and indicates the sensitivity of back- 
scattered infrared observations. We note that such observations pertain to 
the first 3.4 pm optical depth within the cloud, that is, to the upper part of 
the cloud. 

The small value for the minimum detectable number density becomes a large 
value of minimum detectable mass because these particles are so large. Even 
so, it is still interesting to compare the minimum detectable amounts of 
sulfur for both inhomogeneous models. Evaluating equation (39) at the ±25% 
level with the values of F 2 determined above, we find that the minimum 
detectable mass of sulfur in the inhomogeneous cloud model ranges from 17% of 
the mass of acid for 9 pm particles to 38% for 18 pm particles, as compared to 
7% of the mass of acid if the sulfur is within each drop. 


Penetration Height of 3.4 pm Observations 

Finally we estimate the height above the surface of Venus to which these 
3.4 pm observations' and minimum detectable masses pertain. ¥e define the 
effective extinction cross section or coefficient as 


eff.kg(A) = (1 - cos 0)ksa) + k^d) (40) 


where cos 8 is the anisotropy factor (ref. 6). This definition differs from 
that of Lacis (ref. 34, figure 2) in that we apply the correction for forward 
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scattering (1 - cos 0) only to rather than to the sum of kg + k^ as did 
Lacis. Because both the visible polarization and the 3.4 pm observations per- 
tain primarily to single scattered light, we can write 


~^eff.IR 

"^eff.VIS 



eff.kg(IR) N(h)sec z dh 


eff.ke(VIS) N(h)sec z dh 


(41) 


Assuming that eff .k^ sec* z is independent of height gives 



N(h)dh 


eff ,kg(VIS) 
eff.kg(IR) 



N(h)dh 


(42) 


Here N(h) is the number density at height h, and hjg^ and h^ps the 

deepest heights from which radiation of the wavelength specified has arisen 
when is large. 


The column density integral on the right side of equation (42) was deter- 
mined graphically from the density model of Lacis (ref. 34, fig. 8) with 
hvis ~ 66.5 km. The effective extinction cross sections were calculated at 
0.55 and 3.4 ym using the homogeneous model of H^SOj^ drops at 84-1/2% -concen- 
tration and 300 K temperature. Equation (42) was then evaluated graphically 
and produced a value for hpg^ of 67.5 km. We note that Lacis^ density model 
is sufficiently steep at this height that a factor of two reduction in the 
3.4 ym effective extinction cross section (which would result from use' of 
Lacis' definition) would cause only a 2-km Increase in depth. Thus we con- 
clude that the 3.4 ym observations pertain to about the same height as those 
in the visual, or to just slightly higher. 


SUMMARY AND CONCLUSION 


The ability to detect small contributions to the scattering of sunlight 
by sulfur or other nonabsorbing species is greatly enhanced in the infrared 
absorption band of sulfuric acid where absorption reduces backscattering by 
the numerous homogeneous acid drops to levels more than 40 times below those 
in the visual. 

A planetary reflectivity function, whose shape is quite similar to that 
of the single scattering phase function with strong absorption, has been 
defined and related to other known functions. 
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Two inhomogeneous models containing sulfur have been described and tested. 
The principal new parameters associated with these inhomogeneous models are 
the radius Rq of the sulfur core within an inhomogeneous drop, and the 
effective radius a£ and relative number density F 2 of large sulfur par- 
ticles in the inhomogeneous cloud model. A core-size parameter, nl, indicat- 
ing how many inhomogeneous Mie coefficients are required has been defined for 
the inhomogeneous drop model. Calculations with these models at 3.4 ym show 
that backscattering near phase angles of 30° by the constant core, inhomoge- 
neous drop model has a linear dependence on core volume (Rq^); backscatter 
from the inhomogeneous cloud model has a linear dependence on the total area 
of large sulfur particles through either ( 32 )^ or F 2 . The volume of non- 
absorbing core inside the inhomogeneous drop affects .the overall level of the 
3.4 pm reflectivity function by excluding the sulfuric acid absorber; in con- 
trast, the area of the large, non-absorbing, high-index sulfur particles 
affects the shape of the reflectivity function by adding a nearly isotropic 
reflection to the steep curve of the acid drops. 

Minimum detectable amounts of sulfur are based on the ability to observe 
an enhancement of the 3.4 pm backscatter relative to that calculated from the 
homogeneous model. Observations with ±25% or less total error could distin- 
guish the homogeneous drop reflectivity function from that of inhomogeneous 
drops if Rq > 0.4R. Similarly, ±25% observations could distinguish the 
inhomogeneous cloud reflectivity function of sulfur particles with effective 
radii between 9 and 18 pm if F 2 2 . 22x10“^ or 6x10”^, respectively. These 
minimum distinguishable parameters lead to ratios of the minimum detectable ■ 
mass of sulfur to that of H 2 SO 4 for the inhomogeneous drop model of 7%, and 
for the inhomogeneous cloud model of 17% or 38%, depending on particle size. 

These calculations lead us to conclude that moderately accurate 3.4 pm 
observations of the intermediate angle backscatter from Venus are capable of 
detecting quite small amounts of elemental sulfur at heights above the plane- 
tary surface comparable to those of the visual observations. 
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Figure 1.- Reflectivity with strong absorption at a visual wavelength (Hansen 
and Travis's hypothetical case no. 4, ref. 6); P"(<{>) is the single scatter 

ing phase function plotted on an arbitrary scale; R((fi) and p((fi) are the 
reflectivity function and geometric albedo plotted on the same scale. 



Figure 2.- Observation vs theory at 1.0635 iim. The smooth curve is 0.899 
times the calculated reflectivity function; observed data from Irvine 
et al. , (ref. 5) . 
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Figure 4.- Variation of tlie imaginary index in the homogeneous drop model. 
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Figure 5.- Variation of the effective 
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Figure 6.- Variation of the effective variance in the homogeneous drop model. 
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Figure 10.— 'Variation, of real index of- core of inhomogeneous drop model 

slant core calculation. 
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Figure 11.- Comparison of four different ways to vary the core radius of the 
inhomogeneous drop over the size distribution^ 
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Figure 12.- Variation of thickness of constant coating, inhomogeneous drop 

model. 
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Figure 13.- Variation of size of constant core of inhomogeneous drop model. 
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Figure 14.- Variation of size of scaled core of inhomogeneous drop model 
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Figure 15.- Variation of real index of sulfur particles of inhomogeneous cloud 

model. 


ABSOLUTE R {3.4, <f>) 









Figure 17.- Variation of FTB of inhomogeneous cloud model. 
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Figure 18.- Variation of- b 2 of inhomogeneous cloud model. 
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Figure 19.- Variation of effective radius of sulfur particles (a 2 ) of inhoitio 

geneous cloud model. 



Figure 20.- Variation of the fraction particles in the sulfur peak 

of the. ^inhomogeneous cloud model. 
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Figure 21.- Dependence of the backscattered reflectivity function R(3.4,(()g) 
upon sulfur core volume in the inhomogeneous drop model; relative core 
sizes (Rq/R) are indicated by arrows along the abscissa. 
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Figure 23.- Comparison of homogeneous and inhomogeneous models. 
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Figure 25.- Detectability limits of inhomogeneous cloud model for ~ 18 yiii 
and FTB = 3; minimum distinguishable values of F 2 ate indicated by arrows. 
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